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A hydro dynamic description is used to study the zero- 
temperature properties of a trapped spinor Bose-Einstein con- 
densate in the presence of a uniform magnetic field. We show 
that, in the case of antiferromagnetic spin-spin interaction, 
the polar and ferromagnetic configurations of the ground state 
can coexist in the trap. These two phases are spatially seg- 
regated in such a way that the polar state occupies the inner 
part while the ferromagnetic state occupies the outer part 
of the atomic cloud. We also derive a set of coupled hydro- 
dynamic equations for the number density and spin density 
excitations of the system. It is shown that these equations can 
be analytically solved for the system in an isotropic harmonic 
trap and a constant magnetic field. Remarkably, the related 
low lying excitation spectra are completely determined by the 
solutions in the region occupied by the polar state. We find 
that, within the Thomas- Fermi approximation, the presence 
of a constant magnetic field does not change the excitation 
spectra which still possess the similar form of that obtained 
by Stringari. 



The recent realization of Bose-Einstein condensate 
(BEG) in trapped atomic gases has inspired enormous 
interest in the theoretical and experimental studies of in- 
homogeneous interacting Bose gases . More recently, 
Stamper-Kurn et al. have successfully produced a BEG 
of ^■^Na atoms in an optical dipole trap Q| . In their ex- 
periment, a new kind of Bose-condensed system, namely, 
the spinor BEG which is characterized by the three hy- 
perfine spin states |F=l,mi?=:±l,0), has been real- 
ized. A remarkable feature of such a system is that, 
unlike the magnetically trapped BEG, the spin- flip col- 
lisions between atoms allow population exchange among 
hyperfine states without causing any trap loss. This has 
opened up the possibilities to explore the Bose-condensed 
systems having internal degrees of freedom in which the 
C/(l) gauge symmetry as well as the rotational 5*0(3) 
symmetry in spin space are both involved. 

In the spinor BEG, the spin-spin interaction between 
atoms can be approximated by the spin exchange term 
S-S [H where S —Sx^x^ SyGy-'r S z^z is the spin operator. 
In such circumstances, the sign of the spin-exchange cou- 
pling strength may play a crucial role in determining the 
ground state configurations. For example, in the cases 
of negative (ferromagnetic) coupling strength, atoms of 



same spin orientation in the condensate tend to attract 
each other against the repulsive binary scatterings which 
give rise to the density-density interaction. The situation 
is equally complicated for the cases of positive (antifer- 
romagnetic) coupling strength in which atoms of same 
spin orientation tend to repel each other. The competi- 
tion between the density-density and spin-spin interac- 
tions evidently lead to an intriguing scenario of the spin 
dynamics of the spinor BEG that is revealed by certain 
complex ground state structures [^[-^ . 

In the previous theoretical studies, the zero- 
temperature properties of the spinor BEG have been con- 
sidered exclusively either in a trap [^|J^] or in a uniform 
magnetic field [^,1). It is of great theoretical interest to 
ask that how these physical properties are affected when 
the confining trap and the external magnetic field are 
both present. In this paper, we use the hydrodynamic 
treatment to probe into this issue with a special empha- 
sis on the case of antiferromagnetic spin-spin interaction. 
We show that part of our results are related to the for- 
mation of spin domains in the ground state of the spinor 
BEG which has been reported very recently by Stenger 
et al. @. 

We begin by considering the generalized Gross- 
Pitaeviski Hamiltonian 




which describes the classical part of the spinor BEG [g,|6| . 
Here V (r) is the trapping potential and Jl is the Lar- 
mor frequency which is proportional to the magnetic 
field. For simplicity, we shall assume that the magnetic 
field is applied along the z-axis, i.e., Jl — rJ^e^. The 
coupling constants g„ and characterizing the density- 
density and spin-spin interaction, respectively, are given 
by 5n — A'Kan/'rn and gs — A'Kas/m , with a„ and Os 
being the corresponding s-wave scattering lengths. Un- 
like that for the scalar BEG, the order parameter $ for 
the present system is vector-like whose components are 
represented by three classical fields, (pQ and 0_, cor- 
responding to the condensates in the hyperfine states 
\F — 1, mp = 1, 0, —1), respectively. In the presence of a 
magnetic field, there are two classes of symmetry trans- 
formations admitted by the system, namely, the U{1) 



1 



phase symmetry and the SO(2) rotational symmetry 
about the direction of the field. It can be shown that the 
underlying symmetry group of Eq.(|l]) is U {1)^U {1) 
such that the Hamiltonian is invariant under the trans- 
formation: 



-y 



(2) 



Here the superscript T stands for transpose and the pa- 
rameterized phase angles satisfy the relation 9^ + 9^ = 
29o. 

We first examine the ground state configuration for 
the case of ferromagnetic coupling (173 < 0), which is in 
essence identical with that of a scalar BEC. In order to 
minimize the ground state energy, we apply the inequal- 
ity for spin-one systems, i.e., |$^^S4>| < |4>|^ = p(r). Ob- 
viously, the contribution of all the spin-dependent parts 
in Eq.d), 



it is easy to check that the small fluctuation parts d(j)o 
and 6(p- lead to two independent modes whose spec- 
tra are determined by the one-particle Schrodinger equa- 
tions with effective potentials V (r) + giPo — ^ and 
V (r) -I- {gn — gs) Po + — P, respectively. In the homo- 
geneous case, each mode is known to have a free-particle- 
like spectrum with a gap |^,||. The hydrodynamic(-like) 
mode corresponding to the Bose-condensed component 
(f>+ can be obtained by considering order parameter of 
the form 



*(r,t) = (^Jpiv,t)i 



0,0 



(7) 



where p (r, t) — po {r) + Sp (r, t) . Substituting Eq.(0) into 
Eq.(^ and following the analytical approach given by 
Stringari , it is straightforward to show that the den- 
sity fluctuation Sp satisfies the equation of motion 



is minimized by choosing ^^Sz^ 



(3) 



Now since 



6_ I as we have adopted the spin ma- 



trices in the Imp = 1, 0, —1) representation, it follows im- 
mediately that <j>Q = (p- — and ^ 0. Hence the 
ground state configuration based on the hydrodynamic 

description is given by 4>g = ( •\/po(r)e*^''('"\ 0, Oj , 

where po and V^o/w arc identified as the local conden- 
sate density and superfluid velocity, respectively. As we 
wish to deal with a system having fixed total particle 
number N = J cPrpQ (r) , po sjid 9q must be chosen to 
minimize the free energy 



d'r \ {\7^f + ^po (V(?o)' 
' 2m 2m 



+ [V (r) - - /i] PO + \gipt \ , (4) 



where p is the chemical potential, and 171 = gn + gs- 
Obviously, F can be minimized only when V^o = Oi 
that is, ^0 must be position-independent. In the Thomas- 
Fermi (TF) limit, the kinetic energy term (V^/po) /2to 
is neglected and the minimization of the free energy can 
be achieved by letting 



_ J .91 Mm 
0, 



^L-V (r)] , p + nL>V (r) 
elsewhere 



(5) 



We next consider the collective motion of the order pa- 
rameter. In doing so, we assume that the collective modes 
deviate slightly from the ground state by ^$(r,t), i.e., 
$ (r, t) = $g (r)+(5# (r, t). Substituting $ (r, t) into the 
motion equation of the system 
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+gn\^V^ + gs 



(6) 
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V • [poVSp] 



(8) 



Accordingly, the low lying collective excitation spectrum 
is identical with that of the universal excitation spectrum 
obtained by Stringari ||l0[| . 

Of particular interest is the antiferromagnetic case 
{gs > 0). Clearly, the two terms given by Eq.(^ are 
minimized only when — ^^Sy^ — or, ex- 

plicitly, ^X(l)o + (j)*o(t)- = (f>o(t)+ + (t)-(f>Q = 0. Fur- 
thermore, the minimum is achieved if we impose that 
^^Sz^ = ^L/gs > 0. When (/iiQ 7^ 0, by using 
Eq.(H), we may assume without loss of generality that 
(j>o is real. This assumption implies that 0^ — ~(f>- 
and, as a consequence, $^5'^$ = which apparently 
violates the imposed condition ^^Sz^ = flL/gs- We 
thus conclude that 0o must be identically equal to zero 
and are led to the so-called polar state for the condcn- 

T 



sate 

4"^ 



,0, 



where 



: Po (l ± Pags) /2. Note that in this region 
the condition po (r) > ili / gs is required in order to keep 
> 0. On the other hand, if po (r) < ^l/ gs, the 

condition '^'^Sz^ = ^L/gs > can never be fulfilled 
and we must choose $ in such a way that ^^Sz^ is 
as close to n^/gs as possible. In view of the inequal- 
ity, |$^^S4>| < |$| , we conclude that the ground state 
configuration is described by the ferromagnetic state : 



g (r) = (poe'''+",0,0) . We thus see that when 
gs > 0, the two different ground state conflgurations 



$g ' and $g ' coexist and the phase boundary rf, is de- 
termined by Po (fb) = flL/gs- The free energy of such a 
structure of coexistence can be described by 



(fir 



polar 
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{^l/po9s 



2gs 
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fcrro 

+ / d^r 
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(V (r) - p)po + ^gnPl 



(9) 



where we have set VO+q (r) = V^^_o (r) = which serves 
as the necessary condition for minimizing the free en- 
ergy. We subsequently consider the minimization condi- 
tion, SF/Spo — 0, for both regions. In the TF hmit, the 
kinetic energy term is ignored and hence we are able to 
write down the density profile in both regions: 



Po (r) 



g-'[p-V{r) 



polar 



g-^^ [p + — V (r)] , ferromagnetic. 



(10) 



For simplicity, we assume an isotropic harmonic trap 
Vhar (f") = mujQr'^ /2. Thus, the two important length 
scales of the density profile Eg. ([l0|) , namely, the phase 
boundary and the radius of the atomic cloud are 
given by = lUq^ [2 {fj. - ^Lgn/ gs) /m]^^"^ and R = 
u)q^ [2 {fi + ^II) /m]^^'^ , respectively. In this case, one 
sees that the 0_ component is restricted in the region of 
< r < rfc while the (f)^ component exists ubiquitously 
in the atomic cloud. We note that a similar structure 
has been reported in Refs. in the case that the 

quadratic Zeeman energy is almost canceled. The chem- 
ical potential can be obtained by fixing the total particle 
number with the density profile led by Vhar (f) ■ 



N = 



3/2 



15gi yti^Q 

+ (p + 



9s , 
— 

9n 



5/2 



(11) 



This equation provides a complicated mathematical re- 
lation between the chemical potential p and the total 
particle number N. However, when the magnetic field 
is absent, we again obtain the well-known 5/2-power-law 
for the trapped scalar BEC: N — ano C^p/i^o)"^^ /15a„, 
where gho = (™wo) is the characteristic length of 
the harmonic trap. 

In order to obtain the low lying excitation spectrum of 
the configuration described in Eq.(|l0|), we let 



* (r,t) 




polar 

ferromagnetic 
(12) 



Here p± (r,<) = po± (r) + Sp± (r,i), where po± (r) de- 
note the equilibrium particle densities for the fields 
in the configuration of Eq.(^0|) and Sp± (r,t) are the cor- 
responding density fluctuations. In the ferromagnetic re- 
gion the respective linearized equations for the fluctua- 
tions 6p+,S(f>o, and 6(j)- are again decoupled and are, in 



fact, exactly the same as those for the case of ferromag- 
netic coupling. In the polar region, the fluctuation d4>o 
also decouples from 5p±. In fact, we have derived the cor- 
responding linearized equation and found that it is not 
analytically solvable in the presence of a parabolic trap- 
ping potential. Here we shall not discuss that equation 
which is beyond the scope of this paper. Nevertheless, it 
would be instructive to point out that in the homogenous 
case, the motion equation of SipQ is solvable and predicts 
the existence of a massive mode with a gap equal to 
|^,||. Again, we follow the steps in deriving Eq.(||) and 
obtain the following linearized coupled equations: 



V • f ^ [giV5p+ + (.g„ - g,) V5pA , 

L TO J 

V-f^[(g„-5s)V<5p++5iV<5p_]}. (13) 

L TO J 



Let 5p — 5p+ + 5p- and 5ps = 5pj^ — 5p- which 
correspond to the fluctuations of number density and 
spin density respectively. If we further assume that 
the time dependence of these two quantities are de- 
scribed by 5p (r, t) = 5p (r) exp {—iojt) and Sps (r, t) — 
Sps (r) exp (— iwt), we then obtain the desired hydrody- 
namic equations for the system: 



Lu'-'dp - -V • <i gn—ySp + —V5p. 

TO TO 

uj^5ps - -V • ^ ^^^y5p + gAv5p, 

gs TO TO 



(14) 



We note that when the magnetic field is absent, Eq.(|l4|) 
reduces to the same result obtained by Ho On the 
other hand, if the trap is removed, the density po becomes 
a constant and it is thus easy to see that the excitation 
frequencies agree with the results given by Ohmi and 
Machida in the long- wavelength limit . 

In order to solve Eq. (^) for an isotropic harmonic po- 
tential, it is useful to assume solutions of the form: 

5p (r) - / (r) Yira [9, 0) , 5ps (r) = (r) li„, (6, 0) , (15) 

where Yim {0, (f>) is the spherical harmonic function. In 
what follows, we merely substitute the density profile for 
Vhar (r) in the polar region into Eq. ( p^ as it will be illus- 
trated later that this density profile alone is sufhcient to 
solve the excitation spectrum of the system. Introducing 
the dimensionless variable x = r / R, the coupled equation 
for / (a;) and fs {x) in the polar region (Q < x < rjj/ R) is 
expressed by the following matrix equation 



A 

dx"^ 



dx^ X dx 



2d l{l + V) 



Ax 



dx 



f 

f + Bf 



(16) 



0, 



where the matrices are given by 



3 



1 



,B = 



e 

e//3 



,f = 



fix) 

fs{x) 



with the parameters defined by e = 2uP'/u^, x\ = 
2ijOL/mujQR^ and j3 = gs/g-a- Equation ( [iq ) can be 
solved by using the series expansions for the radial func- 

oo 

tions f [x) and fs{x), i.e., f = Ufex'''+' where — 

k=0 

{ak,bk)^ . It is straightforward to show that for k > —1, 
we have, with the convention u_i = 0, the recursion re- 
lation 



Au 



fe+2 



B-k"^ -2kl-ik- 21 
(fc + 2) {k + 21 + 2,) 



(17) 



An immediate consequence of Eq. (^^ is that Ufe vanishes 
for all odd fc. To solve Eq. ( [l7| ) for even k we first work out 
the eigenvalues A± of the matrix A which are given by 
A± = l:kx\/(3. Let u± be the corresponding eigenvectors 
of A± and thus for all fc we have u^. 
Accordingly, Eq. (p^) becomes 



^+4+2"+ + ^- 
B - fc2 - 2kl 



3fc 



21 



(fc + 2) (fc + 2/ + 3) 



(18) 



When /3 7^ 1 ((?„ 7^ g^), it is obvious that u± can not be 
eigenvectors of B. A close examination on the behavior 
of the series in the large order expansion shows that the 
series must terminate at some fc, otherwise it will diverge 
at the phase boundary r = . We therefore conclude 
that for some fc = 2n > 0, the right hand side of Eq.(|l^) 
vanishes and thus the solutions to the eigenvalue prob- 
lems are 



presence of a constant magnetic field does not change the 
frequencies of the hydrodynamic modes but indeed mixes 
up the number density oscillation 5p and the spin den- 
sity oscillation 5ps ■ It is also interesting to point out that 
when gn = gs, we have two independent hydrodynamic 
equation of the Stringari-type which have the same spec- 
tra. 

In conclusion, we have investigated the ground state 
structure of a trapped spinor BEC in the presence of a 
uniform magnetic field. It is found that the antiferromag- 
netic spin-spin interaction leads to the coexistence of the 
polar and ferromagnetic phases in the ground state. For 
the isotropic harmonic trap the density profile is obtained 
in the TF limit. Moreover, the linearized equations for 
the collective modes are derived. We find that the oscil- 
lation frequencies of the hydrodynamic modes are com- 
pletely determined by the linearized equations in the po- 
lar region. Within the TF approximation, the spectra of 
the hydrodynamic modes are shown to be independent 
of the strength of the magnetic field and possess a form 
of Stringari-type. 

Our studies provide a theoretical treatment for the 
zero-temperature properties of the trapped spinor BEC. 
These analyses are, however, limited by the TF approxi- 
mation and the use of an isotropic trap. For all practical 
purposes, it is necessary to extend the present calcula- 
tions to include the effects of the kinetic energy pressure 
as well as the anisotropic trap. 

This work is supported by the National Science Coun- 
cil, Taiwan under Grant NSC-88-2112-M-018-004. 



(a2n, 0)^ , e = + inl + 6n + 21; 



= P (4n2 + Anl + 6n + 21) 



(19) 



In either case, the sequence of the vectors can be de- 
termined by the use of Eq. (^7|) and thus the correspond- 
ing eigenfunctions f can be obtained. For (3=1, the 
situation is simpler: Eq.dT^) becomes two independent 
equations for Sp± . A similar analysis shows that Eq. ( [l9| ) 
is again valid. 

The first mode in Eq.(|l9|) has the spectrum 



,1/2 



■'nl 



loq {2v? + 2nl + 3n + /] ' which possesses exactly the 
same form of the universal spectrum obtained by 
Stringari |lQ| . What is noticeable is the spectrum of the 



1/2 



second mode o;^; = [(2ri^ + 2nl + 3n + /) ^s/ffn] 
which now depends explicitly on the two-body interac- 
tion although it has the same n, I dependence as that of 
^ni ■ When the magnetic field is absent the two modes are 
merely owing to the number density and the spin density 
fluctuations, respectively. It should be noted that both 
two sets of dispersion relations given by Eq.(|l9|) do not 
depend on the magnetic field at all. In other words, the 
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